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Abstract 

We apply the Lie Theory of continuous groups to investigate the sym¬ 
metries of the Hu-Paz-Zhang Master Equation which arises in the 
modelling of the interaction of a harmonic oscillator with a linear 
passive heat bath of oscillators. We examine the case in which the 
parameters of the equation are constant. 
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1 lutroductiou 

Hu, Paz and Zhang derived an (1-1-2) evolution equation with the use of 
path integral methods pen] that inherits the properties of quantum Brow¬ 
nian motion for a harmonic oscillator interacting with a linear passive heat 
bath of oscillators. This equation is widely known as the Hu-Paz-Zhang 
Master Equation. Halliwell and Yu illustrated an alternate and elementary 
derivation of this equation by tracing the evolution equation for the Wigner 
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Function [S]. Ford and O’Connell constructed their results on the basis of 
the quantum Langevin Equation PE]. This equation is given by 

Zt = ~~^y + + 2T{t){xZ)^ + amT{t)h{t)Z^^ + aT(t)f{t)Z^y, 

( 1 . 1 ) 

where Z{t,x,y), a and m are real constants and the time-dependent coeffi¬ 
cients are given by explicit expressions PP. 

The integration of fll.ip is a nontrivial and cumbersome task. Previous meth¬ 
ods that were implemented to determine a solution to fll.ip are either numer¬ 
ical or use a Green’s Function. Neither approach is completely satisfactory. 
The study and analysis of differential equations through the realm of group 
theory is associated with the great mathematician Sophus Lie [T5| and is a 
standard but quite effective technique. The advantage of the Lie theory of 
continuous groups PEaEo] is that it allows for greater generality and, even 
if eventually recourse must be made to numerical methods, this is at the hnal 
stage of the calculation and not ab initio. Such analysis was illustrated in 

PI- 


2 Symmetry analysis 

We analyse the equation given by 

Ru{t, X, y) — xuy + RxUx + ySu^ + Vu^y + Wu^x — «* = 0 (2.1) 

in which there has been a replacement of time-dependent parameters with 
constant parameters and some renaming and rescaling. 

We determine the Lie point symmetries admitted by (12.1 p using the add¬ 
on package SYM PEI SIP associated with the symbolic manipulator Math- 
ematica. 

We take the Lie point symmetry to be of the form 
A = ^^(f, x, y, u)dt + C^(f, x, y, u)dx + x, y, u)dy + r]{t, x, y, u)du (2.2) 
and solve the equation 

{Ru{t, X, y) — xUy -f RxUx -l- ySux -I- Vuxy + Wuxx — Ut) = 0, (2.3) 

where A^^l is the standard second extension of A. 
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The analysis of the determining eqnations is somewhat complicated and 
we simply snmmarise the resnlts. As fl2.ip is linear, it has an inhnite nnmber 
of solntion symmetries which we do not take into consideration. 

The symmetries are 


Ai 

= dt + udu, 


A2 

= udu, 


A3 

= Ai(t)dx + A2{t)dy, 


A4 

= Bi{t)dx + B2{t)dy, 


A5 

= C{t){Cldx + C2dy 

+ {C3X + C^y) udu) 

Ae 

= E{t){Eidx + E2dy 

+ {E3X + E^y) udu) 


where Ai{t), A 2 (t), B 2 {t), C{t), (Ci-Q), E{t), (Ei-Ei) are given by 


Ai{t) 

A2{t) 


'I {R + - 45) exp -^(^R+VR^- 4 :S)t 


exp 


-i (j? + 


B,(() 


i(-B+VB2-4S)exp - (-B + VR? - 4S) ( 


exp 


^ [-R + VR^ - 45 ) t 


C{t) = 

Cl = 
C 2 = 
Cs = 

c, = 


exp 


2{RV + W) 

R - VR'^ - 4SW, 

2 {RV + W), 

R (^-R + VR^ - 45 ) , 
-2RS 


[-R + VR^ - 4A) t 


and 


m = 


2{RV + W) 


exp 


(i? + VR^ - 45 ) t 
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El = R + - ^SW, 

E 2 = 2{RV + W), 

Es = R (-/? - VR^ - 4^) , 
E^ = -2RS. 


Ai represents the autonomy of fl2.ip and A 2 its homogenuity. Neither are 
very suitable for reduction of order. The structures of (As-Ag) indicate that 
they are what has been termed the solution symmetries due to their close 
connection to the solution symmetries of the corresponding Euler-Lagrange 
Equation [H]. In addition we have the relation R"^ — AS > 0. 

In a normal course of events we would replace Ai by simply dt- However, 
it is instructive to look at the reduction using 2dt + Rudu- The transformation 
is 


u{t, X, y) = exp 


' Rt 

.T 


4x,y) , 


(2.5) 


where z satisfies the second-order equation 


Rz xZy RxZx ~\r Sy^x ^xy hh Zxx O' 


( 2 . 6 ) 


3 Reduction of order to an (1+1) equation 


A 3 : The characteristics are 

u,t and r = (^R + VR'^ — dS") y + 2x (3.1) 

and the reduced equation is 

Rz + -{R- - 4S) TZr + 2{y(R+ \/R? - 4S) + 2W) z„ - Zt= a. 

(3.2) 

The symmetries of fl3.2p have the form 


S 


a{t)dt + 


h{t) + -a'{t)r 


dr + z{f{t) + E{r,a,b))dz, 


(3.3) 


excluding the solution symmetries, r)dz, where a, b and / satisfy a third- 
order, a second-order and a hrst-order equation respectively. Thus the re¬ 
duced equation is related to the classical Heat Equation by a point transfor¬ 
mation. 
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A 4 : The characteristics are 


u,t and r = (—-R + V— 45^ y — 2x (3.4) 

and the rednced equation is 

Rz + -[r+ VR^ - 4S) rzr + 2 ( 1 / (fl - \/fl2 - 45) + 2W) z„ - z, = 0. 

(3.5) 

Consequently this also is an (1 + 1 ) equation of maximal symmetry. 

A 5 : The characteristics are 


t,r = KiWy—x and 


u = z exp 


{K^Wy 


x)RKiy 


where 


l-R{K^W + K.,)y^ 

(3.6) 


R-^fW^is , s 

K\ = -;-— and Ko = -r, 

2{RV+ W) {RV+ wy 

and, the reduced equation is 

R (r^ [R - VR‘^ - 4^) +V(^R+ - 4^) + 2W) z + r (v (r^ 

+ RRR^ - 4^) +w{?>R- RR^-AS)) Zr + (vw [r + RR^ - 4^) 

+ 2W^) Zrr - 2 {RV + W)zt = 0. (3.7) 


As in the case of fl3.2p we hnd that the reduced equation has a maximal 
number of Lie point symmetries for an (1 + 1) evolution equation and so is 
a variant of the classical Heat Equation. 

Ag: The characteristics are 


f, r = K^Wy—x and 


u = z exp 


{K^Wy 


x)RK^y 


where 


\r(KIW + K,)y^ 

(3.8) 


R + y/R^ - AS 
2 {RV + W) 


and K 2 is given as above. 
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The reduced equation is 


R (r^ [R + \/i?2 - 45 ) + v[R- VR^ - 4^) + 2iy) z + r (v (r^ 

- RVR^ - 45 ) +w(3R+ VR^-^S)) Zr + (vw (/? - VR‘^ - 45 ) 

+ 2W^) Zrr - 2 {RV + W)zt = 0 (3.9) 

and thus consequently also an (1 + 1 ) equation of maximal symmetry. 

4 Discussion 

We have seen that in all four cases reductions with the solution symmetries 
of ( 12.111 leads to an (1 + 1 ) evolution equation of maximal symmetry and so 
is related to the classical Heat Equation by a point transformation. Inter¬ 
estingly in the reduction from the (1 -|- 2 ) equation to the (1 -|- 1 ) equation 
the number of Lie point symmetries remains the same even though the al¬ 
gebra is different. The 5 -|- 1 -|- 00 symmetries of (12.ip have the algebraic 
structure {Hi W 5 } ©^ ooHi in which the inhnite-dimensional subalgebra 
comprises the solution symmetries of the linear (1 + 2 ) equation ( 12 .ip . the 
one-dimensional abelian subalgebra is Ai and IT 5 is the hve-dimensional 
Weyl-Heisenberg algebra. 

In the case of each of the four solution symmetries, A^-Aq, the reduced (1 + 1) 
evolution equation possesses 3 + 3 + 00 Lie point symmetries with the alge¬ 
braic structure {s/(2,3fj) ©^ W 3 } ©^ cxoHi in which W 3 is the three-element 
Weyl-Heisenberg algebra. Usually one sees 5 + 1 + cxd and the 5 relates to the 
Noether symmetries of the corresponding classical Lagrangian. Here we have 
preferred to write the former structure to emphasis that the hnite algebra 
comprises two subalgebras each of three elements. 

Note that in the above we are using the Mubarakzyanov Classihcation Scheme 
HSl EZl ESI OS] but with some of the algebras listed in the scheme replaced 
by their common physical names thus W 3 is ^. 3^3 and sl{ 2 ,R) is ^. 3 ^ 8 , where 
the latter are the names in the Mubarakzyanov Classihcation Scheme. 

It is interesting to note that the same algebraic structures are to be found 
in the algebraic analysis | 22 ] of the one and two-factor models introduced by 
Schwartz [ 2 l] in his modelling of the pricing of commodities. Consequently 
we have the possibility to transform from the physical situation of a thermal 
bath of oscillators to the commodities market. 
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5 Conclusion 


In this paper we have concentrated upon the autonomous form of the Hu-Paz- 
Zhang Master Equation. Naturally a complete analysis requires examination 
of the nonautonomous form and we expect that to be reported shortly. In 
fact there is a number of problems with explicit dependence upon time which 
are candidates for investigation. These are equations containing parameters 
which are assumed to be constant for ease of analysis, but which, if one 
seeks a closer union with reality, should really be functions of time. A recent 
example of successful analysis is that of the Black-Scholes Equation with 
arbitrary time-dependent coefficients. The algebra was unchanged [23]. One 
could speculate that this would be the case for the Hu-Paz-Zhang Master 
Equation. 
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